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Abstract 

We introduce the Hahn quantum variational calculus. Necessary and sufficient optimality 
conditions for the basic, isoperimetric, and Hahn quantum Lagrange problems, are studied. 
We also show the validity of Leitmann's direct method for the Hahn quantum variational 
calculus, and give explicit solutions to some concrete problems. To illustrate the results, we 
provide several examples and discuss a quantum version of the well known Ramsey model of 
economics. 
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1 Introduction 

Quantum difference operators are receiving an increase of interest due to their applications — 
see, e.g., [IH^. Roughly speaking, a quantum calculus substitute the classical derivative by a 
difference operator, which allows to deal with sets of nondifferentiable functions. In [6], Hahn 
introduced the quantum difference operator Dq^i^, where q g]0, 1[ and uj > are fixed. The Hahn 
operator unifies (in the limit) the two most well known and used quantum difference operators: 
the Jackson g-difference derivative Dq, where q £]0, 1[ (cf. [HIIllHj); and the forward difference A^j, 
where w > (cf. The Hahn difference operator is a successful tool for constructing families 

of orthogonal polynomials and investigating some approximation problems — see, e.g., p^HTB] . 
However, only in 2009 the construction of a proper inverse of Dq^^ and the associated integral 
calculus was given [TflfTS]. 

In this work we introduce the variational Hahn calculus. More precisely, we formulate problems 
of the calculus of variations using Hahn's difference operator and the Jackson-Norlund's integral. 
We discuss the fundamental concepts of a variational calculus, such as the Euler-Lagrange equa- 
tions for the basic and isoperimetric problems, as well as Lagrange and optimal control problems. 
As particular cases we obtain the classical discrete-time calculus of variations [EJ Chap. 8], the 
variational g-calculus [SIE], and the calculus of variations applied to Norlund's sum [20] . 

The plan of the paper is as follows. In Sect. [2] some basic formulas of Hahn's difference 
operator and the associated Jackson-Norlund integral calculus are briefly reviewed. Our results are 
formulated and proved in Sect. [3] Main results of the paper include necessary optimality conditions 
for the basic problem of calculus of variations (Theorem 13.21 in Sect. 13. 1|) and the isoperimetric 
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problem (Theorems 13.31 and 13.41 in Sect. 13. 2[) . as well as a sufficient optimality condition for the 
basic problem (Theorem [23] in Sect. 13.31) . In Sect. 13.41 we show that the direct method introduced 
by Leitmann in the sixties of the XX century [21] . can also be applied with success to quantum 
variational problems via Hahn's difference operator and Jackson-Norlund's integral. Leitmann's 
method is a venerable forty years old direct method that has shown through the times to be an 
universal and useful method in several different contexts — see, e.g., j22l - [34] . Sect. 13.51 provides 
concrete examples of application of our results. Finally, in Sect. 13.61 we apply the developed Hahn 
variational calculus to obtain a quantum version of the well known Ramsey model and we finish 
the paper with conclusions in Sect. S) 



2 Preliminaries 

Let q g]0, 1[ and lu e]0, +oo[ be given. Define j^- Throughout all the paper we assume / 

to be an interval of R containing loq. 

Definition 2.1 (Hahn's difference operator). Let /:/—;■ M. The Hahn difference operator is 
defined by 

provided that / is differcntiable at wo (we are using /'(<) to denote the Frechet derivative). In 
this case, we call Dq_^f the g, w-derivative of / and say that / is g, cj-differentiable on /. 

Example 2.1. Let q g]0, 1[, o; = 0, and 



t^, ifteQ 

if<eM\( 



In this case = 0, and / is g, w-differentiable on the entire real line. However, / is Frechet 
differentiable in zero only. 

Example 2.2. Let q = uj = 1/2. In this case wq = 1. It is easy to see that / : [—1, 1] — > R given by 

' ~t, if t G (-1,0) U (0,1] 
f{t) = {{), ift--l 
1, ift = 

is not a continuous function but is g, w-differentiable in [—1, 1] with 

f-l, if t e (-1,0)U(0,1] 

[-3, ifi = 0. 

Note that, 

limZ?,,„/(t) = A„/(t), \\YaDq,^f{t) = Dqf{t), and lim D,.^f{t) ^ f'{t), 

where 

qt — t 

is the Jackson q-difference derivative [S1|T1[5], and 

(t + 0J) — t 

is the forward difference IQUTTI. 
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Remark 2.1. Let rjq^i^it) :— qt + u, t £ I. Since q g]0,1[ and uj g]0,+cxd[, then 7]q^u}{t) is a 
contraction, rjq^^^I) C /, riq^u{t) <t ior t> wq, Vq,i^it) > t foi" ^ < ^o, and r]q^^{ujo) = t^o- 
The Hahn difference operator has the fohowing properties: 

Theorem 2.1 ( [T71[T5] ). (a) Let f be q,ijj-dijjerentiable on I and Dq^^f = on I. Then f is 
a constant. Conversely, Dq^^c = for any constant c. 

(h) Let f, g he q, io-differentiable at t E L . Then, 

ft) Dq^M + gm = Dq^^fit) + Dq^^g{t), 

(it) Dq^Mam - Dq^M{t))9{t) + f{qt + io)Dq^^g{t), 
ftU) Dq,^ (f) {t) = P™"'*^^^ 9{t)g{qt + c.) / 0. 

(c) f{qt + uj) = fit) + {{qt + t)Dq^^f{t), t e /. 

Example 2.3 ( pTlfTS] ). Let a, e M. We have 

Dq,u.{at + by = aY,Hqt + ^) + bf{at + 6)"-'=-\ (1) 

fe=0 

for n G N and t ^ ojq. 

Following [T7l[18l|35] , we define the inverse of the operator Dq^^ : 

Definition 2.2. Let J be a closed interval of M such that wq, a, 6 G /. For / : / ^ M we define the 
q, oj-integral of / from a to 6 by 



f{t)dq^ut := / f{t)dq^u.t - / f{t)dq,^t, (2) 



where 



f{t)dq^^t [xil -q)-oj)Y, Q'fi^q'' + [k]g,^), ^ e /, (3) 

with [fc]g,tj := '^^iZ'^q ^ for fc e No = N U {0}, provided that the series converges at a; = a and 
X = h. In this case, / is called g, w-integrable on [a.b]. We say that / is q, oj-integrable over / iff 
it is q, w-integrable over [a, 6], for all a, e /. 

Note that, in the integral formulas ^ and ([3]), when w J, we obtain the Jackson g-integral 



f{t)dqt^ / f{t)dqt- / f{t)dqt 



where 



/ fit)dqt = x{l-q)Y,Q''fix<l 

fe=o 



(see, e.g., |36]); while if q t 1 we obtain the Norlund sum 



b ph pa 

f{t)A^t= / /(i)A^t- / f{t)A^t, 



where 

+ 00 



k=o 

(see, e.g., pTIll37ll38j '). This is why the integral defined by © and ^ is called the Jackson-Norlund 
integral. 

If / : / -> K is continuous at wq, then / is g, w-integrable over / (see [l7l[T8] for the proof). 
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Theorem 2.2 ( pnil8p . Assume f : I ^M. be continuous at ljq. Define 

F{x) := / f{t)dg,^t. 



Then F is continuous at loq. Furthermore, Dq_ujF{x) exists for every x £ I and Dq_u]F{x) = f{x). 
Conversely, 



Dq,^f{t)dq^^t = fib) - f{a) 

J a 

for all a,b £ L 

The q, cj-integral has the foUowing properties. 

Theorem 2.3 ( pTlllSj ). (a) Let f,g : I ^ M. be q,uj-integrable on I , a, b, c £ I and k £ 

Then, 

fi) j:f{t)dq.ut^o, 

(ii) sikf{t)dq^^t^ksif{t)dq^^t, 

/a fi^)dq,ujt = - /h° f{t)dq^ujt, 
ftv) J'^ fit)dq,^t = II f{t)dq.^t + ll fit)dq^^t, 
M /a (/ W + 9it)K,^t = J'^ f{t)dq,^t + 9{t)dq.^t. 

(b) Every Riemann integrable function f on I is q,uj-integrable on I. 

(c) If f,g : I -^M. are continuous at ujq, then 



f{t)Dq,,,g{t)dq,,,t = f{t)g{t) 



t=b 



Dq^^{f{t))g{qt + uj)dq,^t, a,b£L 



Property (c) in Theorem 12.31 is the integration by parts formula for the Jackson-Norlund 
integral, and will be useful in the proof of our Theorem 13.21 



Lemma 2.1 ( pTlfTB] ). Let s £ I, f and g be q,uj-integrable. If \ f{t)\ < g{t) for all t £ {g"s 
Hg,w : n £ No}, s £ I, then 











f f{t)dq^^t 


< f 9{t)dq.ut, 


f f{t)dq^^t 


<-j 






J a 


J a 



for all a,b £ {q"s + [n]g^^ : n £ No}. 

However, it should be noted that the inequality 



f{t)dg,^t < / \fit)\dg,^t, 

J a 

a,b £ I, is not always true. For an example we refer the reader to [171 [18]. 

3 Main Results 

Let a,b £ I with a < b. We define the q, cj-interval by 

[a, 6]g,c. := {g"a + W,,^^ : n e No} U + [n],,^ : n e No} U {ua}. 
By V we denote the set of all real valued functions defined on [a, b]q^u: and continuous at loq. 



(4) 
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Lemma 3.1 (Fundamental Lemma of the Hahn quantum variational calculus). Let f G V. One 

has f{t)h{qt + uj)dq i^t — for all functions h £ with h{a) = h{b) ^ if and only if f{t) = 
for all t e [a, b]q^^. 

Proof. The implication "<^=" is obvious. Let us prove the implication "=^" . Suppose, by contra- 
diction, that f{p) ^ for some p S [a, 

Case I. li p ^ wq, then p = q'^a + [k]q^u oi p — q'^b + [k]q^uj for some k G Nq. Observe that 
a(l — q) — id and 5(1 — q) — uj cannot vanish simultaneously. Therefore, without loss of generality, 
we can assume a(l — q) — uj ^ and p = q^a + [k]q^uj- Define 



h{t) = 

Then, 



/(g'^a + [k]q,^), if t = q^+^a +[k + 1],, 
0, otherwise. 



b 

f{t)h{qt + uj)dq,^t = -(a(l - g) - + [k]g,^)h{q''+\ + [k + l]q,^) / 0, 

which is a contradiction. 

Case II. li p = ujq, then without loss of generality we can assume /{ujq) > 0. We know that 
(see [T71IIH] for more details) 

lim g"a + [n\q^ui = hni q'^b + [n\q^ui = uJo- 

n'l'oo nfoo 

As / is continuous at loq, we have 

lim /(q"a + [n]q,^) = lim f{q^b + [n],,„) = f{uJo). 

nyoo njoo 

Therefore, there exists N (zN, such that for all n > the inequalities 

/(q"a+[n],,^) >0 and /(g"6 + [n],,„) > 

hold. If uq ^ a, 6, then we define 

r/(g"6 +[«],,„), if < = g"+ia+[n + l]g,^, for ah n>N 
Ht) = lf{q"a+[n]q^^), if < = + [n + 1],,^, for all n>N 

I 0, otherwise. 



Hence, 



f{t)h{qt + Lo)dq^^t ={b- a)(l -q)Y. ^"/(^"a + Wg,^)/('Z"& + Wg^c) ^ 0, 



which is a contradiction. If cjq = b, then we define 



h{t) 



f{uo), if t = q"+ia+[n + l],,^, for all n > iV 
0, otherwise. 



Hence, 

b 



f{t)h{qt + uj)dq,^t = - / f{t)h{qt + Lu)dq,a = ~{a{l-q)-Lo) ^ q"/(g"a+ [n]g,<,)/(wo) 7^ 0, 



n=N 

which is a contradiction. Similarly, we show the case when luq = a. □ 
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Let E be the linear space of functions y € T) for which the q, w-derivative is bounded on [a, 
and continuous at luq. We equip E with the norm 

\\y\\i= sup \yit)\+ sup \Dq^^y{t)\. 



For s e / we set 



{q"'s + [n]g^^ :neNo}U{a;o}. 



In tire sequel we need one more result. The following definition and lemma are similar to |39| . 

Definition 3.1. Let g : [s]q^^x] — 0, 0[-^ R. We say that g{t, •) is continuous in Oq, uniformly in t, 
iff for every £ > there exists 6 > such that \9 — 6o\ < 6 implies \g{t,6) — 3(^,^0) I < £ for all 
t £ [sjg^tj. Furthermore, we say that g{t, •) is differentiable at Oq, uniformly in t, iff for every e > 
there exists 6 > such that < |(? — ^qI < ^ implies 



'l{t,d)-g{t,0o) 



d2git,9o) 



where d2g — ff for all t £ [s]q,( 



Lemma 3.2. Assume g(t, •) be differentiable at 9o, uniformly in t in [s]q^uj, and that G{9) :~ 
J g{t,9)dq^ujt, for 9 near9o, and J d2g{t,9o)dq,uj exist. Then, G{9) is differentiable at 9q with 
G'\0o)^C^d2g{t,9o)dq^^t. 

Proof. Let e > be arbitrary. Since g{t,-) is differentiable at 9o, uniformly in t, there exists 6 > 0, 
such that, for all t G [s]q^uj, and for < |6' — 0o| < <5, the following inequality holds: 



g{t, 9) - git, 9o) 



-d2git,9o) 



< 



Wo 



Applying Theorem 12.31 and Lemma \TJ\ for < \9 — 0o\ < 5, we have 



G{9) - G(0o) 



G'{9o) 



IL 9{t, 0)dq^^t ~ r g{t, 9o)dq^u.t 



d2git,9o)dq^^t 



-,{t,9) - g{t,9o) 



-d2git,9o) 



dq^U)t 



< 



e e 

dn.ujt 

- Wo 



Wo 



Hence, G(-) is differentiable at 9o and G'{9o) = d2g{t,9o)dq^u:t. 



□ 



3.1 The Hahn Quantum Euler-Lagrange Equation 

We consider the variational problem of finding minima (or maxima) of a functional 



^y]= / f{i,y{qt + i^),Dq.u,y{t))dq^ut 

over all y G E satisfying the boundary conditions 

y{a) = a, y{b) ^ (3, a,(3eR, 



(5) 



(6) 



where / : [a, 6]^^^ X K X M ^ R is a given function. A function y S E is said to be admissi- 
ble iff it satisfies the endpoint conditions ©. Let us denote by d2f and ^3/, respectively, the 
partial derivatives of /(•, •, •) with respect to its second and third argument. In the sequel, we as- 
sume that {u,v) — > f{t,u,v) be a C^(R^,K) function for any t e [a, 6],^^, and f{-,y{-),Dq^^y{-)), 
d2f{',y{'), Dq^c^y{-)), and dzf{-,y{-),Dq i^y{-)) are continuous at wo for any admissible function 
?/(•). Finally, an ft. G E is called an admissible variation provided h{a) — hib) — 0. 
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For an admissible variation h, we define function (f> : ] — e, R by 

0(e) = (t){e\y,h) C[y + eh]. 
The first variation of problem ([S])-® is defined by 

5C[yM :=</.(O;y,/i) = 0'(O). 

Observe that, 

C[y + eh]= I f{t,y{qt + uj) + eh{qt + uj),Dq,^y{t)+eDq^^h{t))dq,^t 

b 

f(t, y{qt + uj)+ eh{qt + uj), Dq,^y{t) + eDq^^h{t))dq^^t 



f{t, y{qt + uj)+ eh{qt + uj),Dq^^y{t) + eDq,^h{t))dq^^t. 



Writing 
and 

we have 
Therefore, 



Cb[y + eh]^ I f(t,y{qt + uj) + eh{qt + uj),Dq^^y{t) + eDq^^h{t))dq^ajt 



Ca[y + eh]^ I f{t,y{qt + uj) + eh{qt + uj),Dq,^y{t) + eDq,^h{t))dq^^t, 



C[y + eh] = Cb[y + sh] ^Ca[y + eh]. 



6C[y,h] = SCb[y,h]~6Ca[y,h]. (7) 
Knowing fT]) , the following lemma is a direct consequence of Lemma 13.21 



Lemma 3.3. Put g{t, e) — f{t, y{qt+uj) + eh{qt+uj), Dq^^y{t) + eDq^^h{t)) for e g] — £, e[. Assume 
that: 

(i) g{t, •) be differentiable atO, uniformly int £ [a]q^i^, andg{t, •) be differentiable atO, uniformly 
in t e [b]q.u,; 

(a) Ca[y + £h] and Cb[y + eh], for e near 0, exist; 

(Hi) J^^d2g{t,0)dq,ujt and ^^^d2g{t,0)dq^ut exist. 
Then, 

5C[y, h]^ [ [d2f{t, y{qt + oj), Dq^^y{t))h{qt + lj) + dafit, y{qt + cj), Dq,^y{t))Dq^^h{t)] dq^^t. 

J a 

In the sequel, we always assume, without mentioning it explicitly, that variational problems 
satisfy the assumptions of Lemma 13.31 

Definition 3.2. An admissible function y is said to be a local minimizer (resp. a local maximizer) 
to problem ([5)) -([6]) iff there exists 5 > 0, such that C[y] < C[y] (resp. C[y] > C[y]) for all admissible 
y with \\y - y||i < 6. 

The following result offers a necessary condition for local extremizer. 

Theorem 3.1 (A necessary optimality condition for problem ([S])-®). Suppose that the optimal 
path to problem ([5])-([6]) exists and is given by y. Then, SC[y,h] — 0. 
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Proof. Without loss of generality, we can assume y to be a local minimizer. Let h be any admissible 
variation and define a function (f> :] — e, e[-^ M by (f)(e) = C[y + eh]. Since y is a local minimizer, 
there exists S > 0, such that C[y] < C[y] for all admissible y with \\y — y\\i < S. Therefore, 
4>{e) = C[y + eh] > C[y] — 4>{0) for all e < p7||7- Hence, has a local minimum at e = 0, and thus 
our assertion follows. □ 

Theorem 3.2 (The Hahn quantum Euler-Lagrange equation for problem ©-(jl])). Suppose that 
the optimal path to problem ([5]) -([6]) exists and is given by y. Then, 

DqMit, y{qt + u:),Dq^^y{t)) = d2f{t, y{qt + w), i^,,^y(t)) (8) 

for all t G [a, h]q_^. 

Proof. Suppose that C has a local extremum at y. Let h be any admissible variation and define a 
function (/):]— e, e[— K by (/)(e) = C[y + eh]. By Theorem 13. 11 a necessary condition for y to be 
an extremizer is given by 

'/''(e)le=o =0^ / [d2f{---)h{qt + ^)+d:,f{---)Dq^^h{t)]dq,^t = Q, (9) 

J a 

where (■••)== (i, y{qt + uj), Dq^^y{t)). Integration by parts (see item (c) in Theorem 12. 3p gives 

dsfi- ■ ■ )Dq^^h{t)dq^^t = 83 f{- ■ ■ )h{t)\lZa - / Dq^Mi- ' ' Mqt + iu)dq,^t. 

Because h{a) — h{b) — 0, the necessary condition © can be written as 

{d2f{- ■■)- Dq^Mi- ■ ■ )) h{qt + Uj)dq^^t 

for all h such that h{a) = h{b) = 0. Thus, by Lemma [3. 11 we have 

d2f{- ■■)- DqMi- • • ) = 

for all t e [a, □ 

Remark 3.1. If the function / under the sign of integration (the Lagrangian) is given by / = 
f{t, yi, . . . , y„, Dq,u,yi, . . . , Dq^ujyn), then the necessary optimality condition is given by n equations 
similar to ([8]), one equation for each variable. 

3.2 The Hahn Quantum Isoperimetric Problem 

Let us consider now the isoperimetric problem, which consists of minimizing or maximizing 

C[y]^ [ f{t,y{qt + Lj),Dq,^y{t))dq^^t (10) 

J a 

over all ?/ € E satisfying the boundary conditions 

y(a) = «, y(5) = /3, (11) 

and the constraint 

J a 

where a, /3, and k are given real numbers. We assume that: 

(i) {u,v) f{t,u,v) and {u,v) g{t,u,v) be C^(R^,R) functions for any t e [a,b]q^uj', 



ny]^ I g{t,y{qt + io),Dq^^y{t))dq,^t = k, (12) 
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(ii) functions /(•, y{-)Dg^^y{-)), d2f{-, y{-), Dq^^y{-)) and 93/(-, y(-), Dg^^y{-)), g{-, y{-),Dq^^y{-)), 
d2g{' ,y{ ), Dq^^y{-)) and d^g{-,y{ ),Dqj^y{-)), be continuous at ujq for any admissible func- 
tion y{-). 

Definition 3.3. An admissible function y is said to be a local minimizer (resp. local maximizer) for 
the isoperimetric problem (ITm - ((T2|) iff there exists (5 > such that C[y] < C[y] (resp. C[y] > C[y]) 
for all admissible y satisfying the boundary conditions (jlip , the isoperimetric constraint p2|) , and 

\\y - y\\i < s. 

Definition 3.4. We say that y is an extremal for IC, iff ^ holds with respect to p^: 

Dq,u:d3g{t, y{qt + uj), Dq,^y{t)) = d2g{t, y{qt + uj), Dq^^y{t)) 

for all t £ [a,b]q^i^. An extremizer (i.e., a local minimizer or a local maximizer) for the problem 
PTI)) - ([T^ . that is not an extremal for IC, is said to be a normal extremizer; otherwise (i.e., if it is 
an extremal for IC), the extremizer is said to be abnormal. 

Theorem 3.3. If y is a normal extremizer for the isoperimetric problem (|10p ~ (|12p . then there 
exists a real X, such that 



Dq^^daFit, y{qt + u), Dq,^y{t)) = d2F{t, y{qt + uj),Dq,^y{t)) 
for all t € [a, where F = f — Xg. 



(13) 



Proof. Consider a variation of y, say y — y + £ihi + £2^2, where hi G E, hi{a) — hi{b) = 0, and £i 
is a sufficiently small parameter, i — 1,2. Here, hi is an arbitrary fixed function and /i2 is a fixed 
function that will be chosen later. Define 



K{ei,e2) ^ IC[y] = / g{t,y{qt + uj), Dq^^y{t))dq.u,t - k. 



We have 



dK 

de2 



(0,0) 



[d2g{- ■ ■ )h2{qt + uo) + dsgi- ■ ■ )Dq^^h2{t)] dq^^t 



where {■•■) — {t, y{qt + uj), Dq^^y[t)). Integration by parts gives 



dK 

de2 



(0,0) 



= / (d2g{---) - Dq^^d3g{---))h2{qt + u!)dq^^t 



since /12(a) 



h2(b) = 0. By LemmaO there exists h2 such that 7^ 0. Since K(0, 0) = 0, 

(0,0) 

by the Dini-Ljusternik implicit function theorem we conclude that there exists a function 62, 
defined in the neighborhood of zero, such that K{ei,e2{si)) = 0, i.e., we may choose a subset of 
variations y satisfying the isoperimetric constraint. 
Let us now consider the real function 

-£(ei, £2) = C[y] = / f{t, y{qt + u), Dq^^y{t))dq^^t. 



By hypothesis, (0,0) is an extremal of L subject to the constraint K = and V-^(0,0) 7^ 0. 
By the Lagrange multiplier rule, there exists some real A, such that V(Z(0,0) — AA'(0, 0)) = 0. 
Having in mind that hi{a) = hi{b) = 0, we can write 



and 



dL 

dei 

dK 
dei 



(0,0) 



= / id2fi---)-Dq^^d3fi---))hliqt+ij)dq,^t 



(0,0) 



= / id2g{- ■■) - Dq^u:d3g{- ■■))hl{qt + U})dq^u,t. 
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Therefore, 



[{d2fi- ■■)- Dg^^dsfi- • • )) - A {d29i- ■■)- Dg^^dsgi- ■ ' ))] hiiqt + uj)d,^^t = 0. (14) 

As HH) holds for any /ii, by Lemma [3. II we have 

52/(- • • ) - ^9-53/(- • • ) - A {d2g{- • • ) - ^9,-535(- • • )) - 0. 
We get 113]) by writing F = f - Xg. □ 

One can easily cover abnormal extremizers within our result by introducing an extra multiplier 
Aq associated with the cost functional. 

Theorem 3.4. If y is an extremizer for the isoperimetric problem (jlOp - (I12[) . then there exist two 
constants Aq and A, not both zero, such that 

Dg,u.d3Fit, y{qt + Lo),D^^^y{t)) = d2F{t, y{qt + iu), D^^^y{t)) (15) 

for all t £ [a, b]q_^, where F = Aq/ — Xg. 

Proof. Following the proof of Theorem l3.31 since (0, 0) is an extremal of L subject to the constraint 
K = the extended Lagrange multiplier rule (see, for instance, [IHl Theorem 4.1.3]) asserts the 
existence of reals Aq and A, not both zero, such that V(Ao£(0,0) — A-ftr(0,0)) = 0. Therefore, 



A 

OEi 



- A — 

(0,0) 



= 

(0,0) 



^ / [Ao (92/(- • • ) - ^g,c.53/(- • • )) - A {d2g{- ■■)- D.^^d^gi- ■ • ))] h^igt + uj)d^^^t = 0. (16) 



Since (|T6|) holds for any hi/\i follows by Lemma 13711 that 

Ao {d2f{- ■■)- Dg^^dsfi- ■ • )) - A {d2g{- ■■)- Dg^^d^gi- • ■ )) = 0. 
The desired condition (fT5|) follows by taking F = Aq/ — Xg. □ 



Remark 3.2. If y is a normal extremizer for the isoperimetric problem ()10p ~ ()12p . then we can 
choose Ao = 1 in Theorem 13.41 and obtain Theorem 13.31 For abnormal extremizers, Theorem 13.41 
holds with Ao = 0. The condition (Ao, X) guarantees that Theorem l3.4l is a nontrivial necessary 
condition. In general we cannot guarantee, a priori, that Ao be different from zero. The interested 
reader about abnormality is referred to the book [41] . 

Suppose now that it is required to find functions yi and y2 for which the functional 

jC-[yi,y2]= f{t,yi{qt + uj),y2{qt + uj),Dq^^yi{t),Dg^^y2{t))dq^^t (17) 

J a 

has an extremum, where the admissible functions satisfy the boundary conditions 

{yi{a),y2{a)) = {yl,y-^) and (2/1(6), (6)) = {y\,yl), (18) 

and the subsidiary nonholonomic condition 

g{t,yi{qt^u),y2{qt + uj),Dq^^yi{t),Dg^^y2{t)) = 0. (19) 

The problem (fT7 |) - (fT9|) can be reduced to the isoperimetric one by transforming (fT9|) into a con- 
straint of the type p^ . For that, we multiply both sides of (fTO)) by an arbitrary function A(i), 
and then take the g, w-integral from a to h. We obtain the new constraint 

JC[yi,y2]^ I X{t)g{t,yi{qt + uj),y2{qt + uj),Dq^^yi{t),Dq^^y2it))dq^^t :^0. (20) 
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Under the conditions of Theorem [231 the solutions {1/1,1/2) of the isoperimetric problem (jl7|) and 
satisfy the Euler-Lagrange equation for the functional 

b 

if - ~X{t)g)dg,^t, (21) 

X{t) = XX{t) for some constant A. Since follows from the solutions of problem ([T7)) - (ITO1) 
satisfy as well the Euler-Lagrange equation for functional (|2T]) . 



3.3 Sufficient Conditions 

In this section, we prove sufficient conditions that ensure the existence of minimum (maximum). 
Similarly to what happens in the classical calculus of variations, some hypotheses of convexity 
(concavity) are in order. 

Definition 3.5. Given a function /, we say that f{t,u,v) is jointly convex (concave) in {u,v), iff 
dif, i = 2,3, exist and are continuous and verify the following condition: 

f{t, u + ui,v + vi) - f{t, u, v) > {<)d2f{t, u, v)ui + dsfit, u, v)vi 

for all {t, u, v),{t, u + ui,v + vi) E [a, b]q,uj x M^. 

Theorem 3.5. Let f{t,u,v) be jointly convex (concave) in {u,v). If y satisfies condition then 
y is a global minimizer (maximizer) to problem ([S])-®. 

Proof. We give the proof for the convex case. Since / is jointly convex in (u, v) for any admissible 
function y + h, we have 

C{y + h)-Ciy) 

b 

[fit, y{qt + + h{qt + io), Dq^^y{t) + Dq^^h{t)) - f{t, y[qt + lo), Dq^^y{tj)] d,,^ 

b 

> I [d2f{t, y{qt + uj),Dg^^y{t))HQt + ^) + dsfit, y{qt + w), Dq^^y{t))Dq^^h{t)] dg,^. 

J a 

We can now proceed analogously to the proof of Theorem I3.2[ As the result we get 
C{y + h) - C{y) > d^f{t, y{qt + uj), mMt)\tl 

+ I {d2f{t, y{qt + uj),y{t)) - Dg^^dsfit, y{qt + w), y{t))) H^t + ^)dg,u,t. 

J a 

Since y satisfy conditions ^ and h{a) ~ h{b) ~ 0, we obtain C{y + h) — C{y) > 0. □ 



3.4 Leitmann's Direct Method 

Leitmann's direct method permits to compute global solutions to some problems that are vari- 
ationally invariant under a family of transformations |2ll[26l[27l|32l[33] . It should be mentioned 
that such invariance transformations are useful not only in connection with Leitmann's method 
but also to apply Noether's theorem [1^143) . Moreover, the invariance transformations are related 
with the notion of Caratheodory equivalence [33J|33] . 

Recently, it has been noticed by the authors that the invariance transformations, that keep the 
Lagrangian invariant, do not depend on the time scale |31j . This is also true for the generalized 
Hahn quantum setting that we are considering in this work: given a Lagrangian / :]R.xMx]R^M, 
the invariance transformations, that keep it invariant up to a gauge term, are exactly the same 
if the Lagrangian / is used to define a Hahn quantum functional ([5]) or a classical functional 
C[y] ~ L{t, y{t) , y' {t))dt of the calculus of variations. Thus, if the quantum problem we want 
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to solve admits an enough rich family of invariance transformations, that keep it invariant up to 
a gauge term, then one need not to solve a Hahn quantum Euler-Lagrange equation to find its 
minimizer: instead, we can try to use Leitmann's direct method. The question of how to find the 
invariance transformations is addressed in [55JI1S]- 

Let / : [a,6]g_i^ x R x R — )■ R. We assume {u,v) f{t,u,v) is a C^(R^,]R) function for any 
t e [a, and f{-,y{-),Dq^^y{-)), d2f{-,y{-), Dq^^y{-)), and d3f{-,y{-),Dg^^y{-)) are continuous 
in ujQ for any admissible function y{-). Consider the integral 

C[y]^ I fit,y{qt + oj),Dq^^y{t))dq^^t. 

Lemma 3.4 (Leitmann's fundamental lemma via Hahn's quantum operator). Let y — z{t,y) be 
a transformation having an unique inverse y = z{t,y) for all t G [a,b]qj^, such that there is a 
one-to-one correspondence 

y(t) ^ m 

for all functions 2/ G E satisfying ([6]) and all functions y G E satisfying 

y^z{a,a), y = z{b,(3). (22) 

// the transformation y = z{t, y) is such that there exists a function G : [a, bjq^^j x R — > R satisfying 
the functional identity 

fit, y{qt + Lo), Dq,^y{t)) - f{t, y{qt + u),Dq^^yit)) = Dq^uGit, y{t)) , (23) 

then if y* yields the extremum of C with y* satisfying (j22p . y* = z{t,y*) yields the extremum of 
C for y* satisfying 1^ . 

Remark 3.3. The functional identity (|23p is exactly the definition of variationally invariance when 
we do not consider transformations of the time variable t (cf. (4) and (5) of [33]). Function G 
that appears in ([23| is sometimes called a gauge term [43]. 

Proof. The proof is similar in spirit to Leitmann's proof [2Tl[26l[27l[30] . Let y G E satisfy ([6]) , and 
define functions y G E through the formula y — z{t, y), t G [a, b]q^ui- Then y G E and satisfies (1^^ . 
Moreover, as a result of ((23]) . it follows that 

pb j*b 

C[y] - C[y] = / /(<, y{qt + uo), Dq,^y{t))dq^^t - / /(i, y{qt + u), Dq^^y{t))dq^^t 

J a J a 

b 

Dq^^Git, y{t))dq.^t = G(&, m) - G{a, y{a)) 

= G{b,z{b, f3)) ~G{a,z{a,a)), 

from which the desired conclusion follows immediately since the right-hand side of the above 
equality is a constant, depending only on the fixed-endpoint conditions (jG]). □ 

Examples 13. 2[ 13.31 and 13.41 in the next section illustrate the applicability of Lemma 13.41 The 
procedure is as follows: (i) we use the computer algebra package described in [45] and available from 
the Maple Application Center at |http: //www.maplesof t . com/applications/view. a spx? SID=4805| 
to find the transformations that keep the problem of the calculus of variations or optimal con- 
trol invariant; (ii) we use such invariance transformations to solve the Hahn quantum variational 
problem by applying Leitmann's fundamental lemma (Lemma 13. 4p . 
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3.5 Illustrative Examples 

We provide some examples in order to illustrate our main results. 

Example 3.1. Let q, lu be fixed real numbers, and / be a closed interval of M such that cjq, 0, 1 G /. 
Consider the problem 

minimize C[y] = ^ (^y{qt + lj) + ^{Dq,^y{t)f^ dq^^t (24) 
subject to the boundary conditions 

y(0)=0, y(l) = l. (25) 

If y is a local minimizer to problem (P^ - (l^5]) . then by Theorem l3.2l it satisfies the Euler-Lagrange 
equation 

Dg,^Dg^^y{t) = 1 (26) 

for all t e {|^[?^]<^ : n G No} U {q'^ +uj[n]q : n G No} U {wo}- By direct substitution it can be verified 
that y{t) = ^q^L^^ ~ gTT^ ^'^ ^ candidate solution to problem 

In next examples we solve quantum variational problems using Leitmann's direct method (see 
Sect, [la- 



Example 3.2. Let g, a;, and a, b (a < b) be fixed real numbers, and / be a closed interval of M 
such that Wo G / and a, 6 G {q"s + [n]q^u : n G No} U {ujq} for some s G /. Let a and /3 be two 
given reals, a ^ 13. We consider the following problem: 



(27) 



.6 

minimize C[y\ = / {{Dq^^y{t)Y + y{qt + w) + tDq^^y{t)) dq^^jt , 

J a 

y{a)^a, y{b) ^ /S . 
We transform problem (1271) into the trivial problem 

minimize C[y] ^ / {Dq^^y{t)Ydq^u,t ■, y{a) = , y{b) ^ , 

J a 

which has solution y = 0. For that we consider the transformation 

y{t) = y{t) +ct + d, c, d G K, 
where constants c and d will be chosen later. According to the above, we have 

Dq,ujy{t) = Dq^^y{t) + c, y{qt + w) = y{qt + u) + c{qt + u>) + d, 

and 

{Dq^^y{t)f + y{qt + uj) + tDq^^yit) 

= {Dq,u,y{t)f + 2cDq^^y{t) +c^ + y{qt + co) + c{qt + co) + d + tDq^^y{t) + ct 
= {Dq^^y{t)f + Dq^^[2cy{t) + ty(t) + ct' + (c^ + d)t\. 

In order to obtain the solution to the original problem, it suffices to chose c and d so that 

ca + d = a, cb + d^ (3. (28) 

Solving the system of equations (|28p we obtain c = and d — . Hence, the global 

minimizer for problem ([27]) is 

a ^ B [3a ^ ba 

y{t) = '^t+'- -. 

a — a — 
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Example 3.3. Let q, uj, and a, b {a < b) be fixed real numbers, and / be a closed interval of R 
such that e / and a,b ^ {^"s + [n]g^^ : n £ No} U {wq} for some s ^ I. Let a and /3 be two 
given reals, a ^ j5. We consider the following problem: 



2 

mmimize C[y] ^ [Dq^^{y{t)g{t))] dq^^t , y{a) = a , y{b) = (3 , 

J a 



(29) 



where g does not vanish on the interval [a, b]q^uj- Observe that y{t) — g ^{t) minimizes C with end 
conditions y{a) = g^^{a) and y{b) = g^^{b). Let y{t) — y{t) +p{t). Then 

[DqMt)9it))f = [Dq.^imgit))? + DqMt)9{t))Dq,^ {m)9it) +p{t)g{t)) . (30) 

Consequently, if p{t) = {At + B)g~^{t), where A and B are constants, then (15(11) is of the form 
(1^51) . since Dq^^{p{t)g{t)) is constant. Thus, the function 

y{t)^{At + C)g-\t) 

with 

A = [ag(a) - Pg{b)] {a - b)-\ C = [aPg{b) - bag{a)] {a - b)-\ 

minimizes (|29p. 

Using the idea of Leitmann, we can also solve quantum optimal control problems defined in 
terms of Hahn's operators. 

Example 3.4. Let w be real numbers on a closed interval / of R such that luq G I and 0, 1 S 
{q"'s + [n\q^i^ : n G No} U {i^o} for some s € I. Consider the global minimum problem 

minimize £[ui,U2]= / {{ui{t))'^ + U2(t))'^) dq^^^t (31) 
Jo 

subject to the control system 

Dq^^yiit) ^ eiip{ui{t)) + ui{t) + U2{t) , Dq^^y2{t) = U2{t) , (32) 

and conditions 

2/i(0) = 0, 2/i(l) = 2, 2/2(0) -0, 2/2(1) = 1, ui{t) ,U2{t) e n = [-1,1] . (33) 

This example is inspired from |33l . It is worth to mention that due to the constraints on the values 
of the controls {ui{t), U2{t) S ri = [—1, 1]), a theory based on necessary optimality conditions to 
solve problem (|3T|) -(p3 l) does not exist at the moment. 

We begin noticing that problem (|3T |) ~(p3 l) is variationally invariant according to [45] under the 
one-parameter family of transformations 

yl^yi+st, 2/2 = 2/2 + si , = ^2 + s {t" and ul = ui) . (34) 

To prove this, we need to show that both the functional integral C and the control system stay 
invariant under the s-parameter transformations ([34]). This is easily seen by direct calculations: 



"2J 



{Ul{t)f + (Umf dq^^t 

1 

Ul{t)^ + {U2{t) + Sf dq^^t 

(35) 

1 

{ui{tf + U2{tf + Dq^^t[sh + 2sy2{t)]) dq^^t 
U2] + + 2s . 
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We remark that and C have the same minimizers: adding a constant + 2s to the functional 
C does not change the minimizer of C. It remains to prove that the control system also remains 
invariant under transformations (j34p : 

Dq.^ (ytit)) = Dg^^ {yi{t) + st) = Dq^^yi + s 

^ exp(u^(t)) + u\{i) + ul(t) , 
Dq,^ ivlit)) = Dg^^ (2/2(0 + St) = Dq^^y2 + s 

= um ■ 

Conditions (1551) and p6l) prove that problem pip - ([55)) is invariant under the s-parameter transfor- 
mations p4l) up to Dq^i^ i^s^t + 2sy2{t)). Using the invariance transformations ([M)) . we generalize 
problem pip - (l33l) to a s-parameter family of problems, s G R, which include the original problem 
for s = 0: ^ 

minimize C''[ui,U2] = / (uKt))"^ + {u2{t))'^dq,^t 
Jo 

subject to the control system 

Dq^^ (ym) = exp{ul{t)) + ulit) + uf,{t) , Dq,^ {y'2it)) = "2 W , 

and conditions 

1/^(0) = 0, yj(l)=2 + s, 2/^(0) =0, y2^(l) = l + s, 

Mi(i) e [-1,1] , ^2(0 e [-1 + 5,1 + s]. 

It is clear that > and that = if uf (i) = u|(i) = 0. The control equations, the boundary 
conditions and the constraints on the values of the controls imply that uf{t) = M|(t) = is 
admissible only if s = —1: yl~^^{t) = t, y2~~^it) = 0. Hence, for s = — 1 the global minimum to 
is and the minimizing trajectory is given by 

Since for any s one has by (|35|) that £[wi, M2] — -C*[uf , — — 2s, we conclude that the global 
minimum for problem £[ui, U2] is 1. Thus, using the inverse functions of the variational symmetries 
(El, 

ui{t)=ul{t), U2{t) ^ ul{t) ~ s , yiit) ^ ylit) - St , y2{t) ^ yl{t) - st , 
and the absolute minimizer for problem pip ~ (l33l) is 

wi(t)=0, U2{t)^l, yi{t)^2t, y2{t)^t. 

3.6 An Application Towards a Quantum Ramsey Model 

As the variables, that are usually considered and observed by the economist, are the outcome of a 
great number of decisions, taken by different operators at different points of time, it seems natural 
to look for new kinds of models which are more flexible and realistic. Hahn's approach allows for 
more complex applications than the discrete or the continuous models. A consumer might have 
income from work at unequal time intervals and/or make expenditures at unequal time intervals. 
Therefore, it is possible to obtain more rigorous and more accurate solutions with the approach 
here proposed. 

We discuss the application of the Hahn quantum variational calculus to the Ramsey model, 
which determines the behavior of saving/consumption as the result of optimal inter-temporal 
choices by individual households [47j . For a complete treatment of the classical Ramsey model we 
refer the reader to i48j . Before writing the quantum model in terms of the Hahn operators we will 
present its discrete and continuous versions. The discrete-time Ramsey model is 

T-l 

max V'(H-p)~*J7 



exp(wi(i)) + ui[t) + U2[t) + s 



U2{t) + s 



(36) 



t+1 



Wt+1 
1 + r' 
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while the continuous Ramsey model is 



max 

W(-) 



-P*U [rW{t) - W'it)] dt, C{t) = rW{t) - W'{t), 



(37) 



where the quantities are defined as 

• W - production function, 

• C - consumption, 

• p - discount rate, 

• r - rate of yield, 

• U - instantaneous utility function. 

One may assume, due to some constraints of economical nature, that the dynamics do not 
depend on the usual derivative or the forward difference operator, but on the Hahn quantum 
difference operator Dq^^. In this condition, one is entitled to assume again that the constraint 
C{t) has the form 



C{t) 



E -r, 



t — UJ 



E 



t-oj 



Wit) 



where E (2, •) is the g, w-exponential function defined by 

00 

E{z,t) ll{l + zq''{t{l-q)-io)) 



k=0 



for z G C. Several nice properties of the g, w-exponential function can be found in |17[|18j . By 
taking the 5, cj-derivative of E ^(^) following is obtained: 



C{t) 



E 



t — UJ 



E\-r,- 



t — UJ 



Dq.u,W[t) 



E[-r,- 



t — UJ 



W{qt + uj) 



(l + r{t-'-^)){l~r{t{l~q)-Uj)) 



The quantum Ramsey model with the Hahn difference operator consists to 



max 

W(-) 



E{~p,t)U 



W{qt + uj) 



r 1 + r U - 



- r 1 - 



(l + r(t-t^)){l-r{t{l-q)-uj)) 



Dq,^W{t) 



subject to the constraint 

C{t) = W{qt + uj) 



r 1 + r U 



t — UJ 



r 1 



1 + rit- 



{1 - r it{l - q) - uj)) 



(38) 



(39) 



The quantum Euler-Lagrange equation is, by Theorem 13. 2[ given by 



(l + r(t-i^))(l-r(t(l-g)-c.)) 



Dq,^[E{-p,t)U'[C{m=0. (40) 



E{-p,t)U' [C{t)] 



Note that for q '[ I and w 4- problem (p8)) -(p9 l) reduces to (|37p . and (|40|) to the classical Ramsey's 
Euler-Lagrange differential equation. 
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4 Conclusion 



In this paper we consider variational problems in the context of the Hahn quantum calculus. Such 
variational problems are defined through the Hahn quantum difference operator and the Jackson- 
Norlund integral. The origin of the Hahn quantum difference operator dates back to a 1949 paper 
of W. Hahn [6] where it was introduced to unify, in a limiting sense, the Jackson g-difference 
derivative and the forward difference. For both of these latter two quantum difference operators, 
variational problems have been studied previously. The forward difference problems were studied 
at least as early as 1937 by T. Fort [3D] and for the g-difference by G. Bangerezako in 2004 [2]. 
In both of these works the authors discuss necessary conditions for optimality and obtain the 
analogue of the classical Euler-Lagrange equation, as well as other classical results. The goal 
of the present paper is to provide extensions of the previous results for the more general Hahn 
quantum difference operator. 

Another related course of study is that of the notion of a time scale. The origins of this idea 
dates back to the late 1980's when S. Hilger introduced this notion in his Ph.D. thesis (directed by 
B. Aulbach) and showed how to unify continuous time and discrete time dynamical systems |49] . 
Since this important result, the literature has exploded with papers and books on time scales in 
which many known results for ordinary differential equations and difference equations have been 
combined and extended [SIllMlllZ]- The classical results of the calculus of variations have been 
extended to times scales by M. Bohner in 2004 [S^. However, the Euler-Lagrange equation here 
obtained is not comparable with that of [32] ■ Indeed, the Hahn quantum calculus is not covered 
by the Hilger time scale theory. This is well explained, for example, in the 2009 Ph.D. thesis of 
Aldwoah ^7] (see also [35]). Here we just note the following: if in Bohner 's paper [39] one chooses 
the time scale to be the g-scale T := : n G Z}, then the expression of the delta-derivative 
coincides with the expression of the Jackson g-difference derivative. However, they are not the 
same. There is an important distinction: the Jackson g-difference derivative is defined in the set 
of real numbers while the time-scale derivative is only defined in a subset T of the real numbers. 
One more difference, between the Hahn calculus we use in this paper and the time scale theory, 
is the following: the delta integral satisfies all the usual properties of the Riemann integral while 
this is not the case with the Jackson-Norlund integral: the inequality ([4|) is not always true for 
the Jackson-Norlund integral. 

The main advantage of our results is that they are able to deal with nondifferentiable functions, 
even discontinuous functions, that are important in physical systems. Quantum derivatives and 
integrals play a leading role in the understanding of complex physical systems. For example, 
in 1992 Nottale introduced the theory of scale-relativity without the hypothesis of space-time 
differentiability [50]. A rigorous mathematical foundation to Nottale's scale-relativity theory is 
nowadays given by means of a quantum calculus [HSj. We remark that results in Bohner's 
paper |39j are not able to deal with such nondifferentiable functions. Variational problems in j39) 
are formulated for functions that are delta-differentiable. It is well known that delta-differentiable 
functions are necessarily continuous. This is not the case in our context: see Example 1 2. 2 1 where 
a discontinuous function is g, w-differentiable in all the real interval [—1, 1]. 

We believe that the obtained results are of interest in Economics. Economists model time as 
continuous or discrete. The kind of "time" (continuous or discrete) to be used in the construction 
of dynamic models is a moot question. Although individual economic decisions are generally made 
at discrete time intervals, it is difficult to believe that they are perfectly synchronized as postulated 
by discrete models. The usual assumption that the economic activity takes place continuously, is a 
convenient abstraction in many applications. In others, such as the ones studied in financial-market 
equilibrium, the assumption of continuous trading corresponds closely to reality. We believe that 
our Hahn's approach helps to bridge the gap between two families of models: continuous and 
discrete. We trust that the field here initiated will prove fruitful for further research. 
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